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We show that the non-standard neutrino interactions can play a role as sub-leading effect on 
the solar neutrino oscillations. We observe that very small flavor universality violations of order 
of O.f-0.2 Gf is sufficient to induce two phenomena: suppression of the i/ e -earth regeneration and 
a shift of the resonance layer in the sun. We obtain these phenomena even in the absence of any 
flavor changing interactions. We discuss their consequences and confront with a global analysis of 
solar+KamLAND results. We conclude that a new compatibility region in the Am 2 x tan 2 8q , 
which we call very low Large Mixing Angle region is found for Am 2 ~ 10" 5 eV 2 and tan 2 = 0.45. 

PACS numbers: 14.16.Pq, 28.50.Hw 



I. INTRODUCTION 

In the last years, the discovery of neutrino oscillation in solar and reactor experiments selected as a more probable 



explanation to the solar neutrino problem the so called Large Mixing Angle (LMA) MSW solution. The SNO 1, 
and the KamLAND ^ experiments confirm and refine the trend of the evidences of neutrino oscillations due the solar 
neutrino observations, as measured by Homestake SAGE GALLEX Q, GNO |jj and Super-Kamiokande |sL flpf . 
As a result, the solar oscillation parameters have pinned down to 6 x 10~ 5 eV 2 < Am 2 < 1 x 10 -4 eV 2 and 
0.3 < tan 2 8q < 0.55 at 2a Several analyzes have arrived to same conclusions 

In a more general context, sub-leading effects can change this picture, which motivate us to investigate the robustness 
of the determination of the solar parameters. In this letter, we assume that non-standard neutrino interactions, which 
we parameterized by two parameters e and e, are present, relaxing the allowed region of the parameters. In the 
presence of non-standard neutrino interactions, we have found that the allowed interval for Am 2 increases, rescuing 
the very low LMA region, Am 2 ~ 1 x 10 -5 eV 2 , and the high part of LMA region, Am 2 ~2x 10~ 4 eV 2 , respectively 
due the suppression of earth matter and due to a Am 2 shift induced by a non-zero e' . 

II. NEUTRINO EVOLUTION: MSW MECHANISM AND NON-STANDARD NEUTRINO 

INTERACTIONS 

We will work in a two generation neutrino scheme, with the contribution from the non-standard neutrino interactions 
(NSNI) QflQ added to the usual MSW Hamiltonian Q. The Hamiltonian in the flavor basis equals 

H = Hmsw + Hnsni 
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where 
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V2G F e f N f { 
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where Nt = N e + 2N n when the NSNI occur with d-quarks, Nf = 2N e + N n when u-quarks are involved and simply 
Nf = N e when we have electrons. The parameters e and e' describe, respectively, the relative strength of the flavor 
changing neutrino interactions and the new flavor diagonal, but non-universal interactions, normalized to G F . 

The NSNI parameters are constrained by non-universal and flavor-changing processes to be e' < 0.7 and e < 

r~i 

10~ 2 20]. Since we concentrate on the regions around LMA neutrino parameters, effects of e parameters will be 
negligibly small and we can, effectively, set this par ameter to zero. We solve numerically the evolution equation, using 
the density profile of the Sun |2| or the earth . 

We now discuss the behaviour of ^-survival probability that will help to understand our results. In the usual MSW 
mechanism, for the solar parameters in the LMA region, we have a resonant behaviour inside the sun given by 



A ™ 2 nn \ ( 2 \ Am 2 cos 2(9 

cos 29 — = — = 
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In the Sun, N n ~ [0.1 — 0.3]iV e , and a positive value of e' can be interpreted as a small negative correction in the 
solar density. For a given Am 2 and 9, the resonance is displaced to the center of the Sun, and as a consequence, less 
neutrinos experience the resonance. As a result, the transition between resonant and non-resonant survival probability 
is displaced to higher values of Am 2 /AE, according to: 

Am 2 Am 2 1 1 „ . Ar . . 

> -—=G F e' f Nf(r) . 

AE AE cos 26> ^2 S 1 

As stated above, effects of a non- vanishing e are much weaker, since for LMA we have Ara 2 /(4i?) ~ 10 -12 eV, 
the same order of \/2G F Nf(r) at the resonance region in the sun. Therefore we should have e ~ 1, much above the 
experimental limit, to have some effect on the survival probability. We neglect e in what follows. 

In the earth N e = N p ~ N n ~ l/ZNd- Therefore, for values of e' ~ 1/3 the matter term in the evolution matrix 
due to NSNI has the same order of magnitude of the ordinary matter term. As a consequence, the regeneration of v e 
is suppressed. 

Both these effects are presented in Fig. ^ Around Am 2 /AE ~ [10 -12 — 10~ n ] eV we can see the displacement in 
Am 2 /AE of the survival probability in the Sun, and around Am 2 /AE ~ [10 -14 — 10 -12 ] eV 2 the suppression in the 
regeneration effect is effective. 



III. SOLAR NEUTRINO AND KAMLAND DATA ANALYSIS 

We use for the solar neutrino analysis the same data set and the same procedure of analysis appearing in Ref. 
Here the main ingredients of the analysis are summarized. 
The data sample consists of 

- 3 total rates: (i) the Ar-production rate, Qat, from Homestake [jj], (ii) the Ge— production rate, Qc e from SAGE 
Q and (hi) the combined Ge— production rate from GALLEX and GNO 



- 44 data points from the zenith-spectra measured by Super-Kamiokande during 1496 days of operation 

- 34 day-night spectral points from SNO plus CC, NC and ES rates from SNO salt-phase Q; 

- 3 fluxes from the SNO salt phase |3| measured by the CC-, NC and ES- reactions. 

Altogether the solar neutrino experiments provide us with 84 data points. All the solar neutrino fluxes are taken 
according to SSM BP2000 Q 

Thus, in our analysis of the solar neutrino data we have three fit parameters: Am 2 , tan 2 9q and e'. 

We define the contribution of the solar neutrino data to x 2 as 

X 2 sun= E (Kih-RLKfiRL-Rix) , (2) 
t,i=l,84 

where we construct the 84 x 84 covariance matrix afj taking in consideration all correlations between uncertainties. 

Following the procedure done in Ref. [2^. the KamLAND data are analyzed through a Poisson statistics, using the 
following x 2 : 

(pjobs 

where the In term is absent when bins with no events are considered (5 last bins). 

The combined analysis of solar + KamLAND data is done just adding the two contributions in x 2 '- 



2 _ 2 1 2 

X XKL ' Xsun 

We minimize the global x 2 with respect to the three parameters Am 2 , tan 2 8 Q and e'. For the KamLAND x\l^ ^ ne 
effect of e' is negligible due to the short distance traveled inside earth, then effectively x\l depends only on Am 2 and 
tan 2 6q. We show our results in the plane Am 2 and tan 2 9q, in Fig. [21 where we minimized away the dependence on 
the e' parameter. The best fit is for d — 0, and the usual MSW mechanism is still the best solution and any non-zero 
value for e' parameter only worsens the fit. 

To have an idea of the bounds put by the present data and the possible accumulation of more statistics on Kam- 
LAND, we plot in the upper panel of Fig. [2Jthe bounds on d parameter. The present bound is shown by a thick solid 



curve, where we have e' d < 0.3 at 2a. Assuming a future exposure of 1 kton-yr for KamLAND experiment and the 
present solar neutrino data, we simulate the KamLAND data as generated by a specific Am 2 and tan 2 8q combination 
located in each one of three islands: very low, normal, and high LMA region. If we assume a point in the normal LMA 
region, we will get after 1 kton-yr for KamLAND, —0.4 < e' d < 0.25 at 2a as showed in the dashed curve. Similar 
plots for the high (very low) LMA region as long-dashed (dotted) curves show that the bounds will be e' d < —0.16 
(0.16 < e' d < 0.35). 

A. NSNI with u-quarks and electrons 

All our results showed were computed for a NSNI with d-quarks. In this section we summarized the main differences 
when you have NSNI with u-quarks and electrons. 

To get NSNI with u-quarks, we need to replace Nd — > N u in the evolution equation for neutrinos. In the production 
region where most of 8 B neutrinos are created, the ratio between the Nd and N u densities is almost constant, as 
can be seen in Figure 3 of Ref. j^jj. Then the conversion probability showed in Fig. ^ for d-quarks, is similar to 
the conversion induced by u-quarks, with the appropriate rescaling of e' u parameter: e' u — e' d * Nd/N u . In the earth 
the differences are minimal. If you compare the allowed region for d-quarks (see Fig. |3 ) and u-quarks then we have 
practically the same figure for u-quarks. 

For NSNI induced by electrons, we can think as a rescaling of the usual matter potential of the MSW mechanism: 
\f2GvN e (r) -> V2G F N e {r) - y/2G F e' e N e (r) = y/2(l - e' e )G F N e (r). The parameter e' e have looser bounds then 
e'd |2j| and values for e' e ~ 1 are still allowed. One could worry that such large values of e' e could cause a strong 
effect on the detection cross section a(ve — ► ve) used to detect solar neutrinos on Super Kamiokande experiment. 
This is not true because the matter potential induced by e' e is proportional only to the vector contribution of the 
non-standard neutrino couplings, and the cross section depends on a combination of the left/right couplings of the 
non-standard neutrino interaction of the electron neutrino. 

When we looked for the allowed regions, we have similar figures as in Fig. [3 and we obtain again the appearance of 
a very low LMA region. We have also shown in Fig. OH (second and third panels) the limits in e' for NSNI induced by 
u-quarks and electrons. For comparison, similar plot was obtained in Ref. [3], bu. 

IV. DISCUSSION OF RESULTS 

In the allowed region showed, we notice two distinguished facts: the appearance of new very low LMA region at 
Am 2 ~ 10~ 5 eV 2 and the old very high LMA, Am 2 ~lx 10 -4 eV 2 , appears at la. 

The analysis of the data, in the absence of d parameter, disfavor the high LMA solution that is only allowed at 
3cr |ll| . Assuming a negative e', the situation changes due the shift of the resonance layer in the sun, as commented 
before. For e' = —0.6 (the lower limit we used in our analysis) this region is now accepted at la C.L.. 



The most interesting phenomenology happens at the island of allowed region in KamLAND analysis around Am 2 ~ 
10~ 5 eV 2 . This region is not allowed in a pure-LMA scenario due to the very high day-night asymmetry that is expected 
for these parameters. But when we include a e' > 0, the NSNI term in Hamiltonian compensates the effect of Earth 
matter interaction. For e' ~ 1/3 we expect a very low regeneration, and for larger values of e' we can have even 
a positive day-night asymmetry. This region is also allowed when we assume other non-standard mechanisms, as 
random matter density fluctuations in the sun, as presented in [2^ . 

V. CONCLUSIONS 

We showed that NSNI will affect the fit in the LMA region of the MSW solution to the solar neutrino anomaly. 
When one takes into account the KamLAND results, positive values of the e' push the allowed region of the neutrino 
parameters Am 2 and tan 2 8q at 95% C.L. from pure MSW low-LMA and high-LMA to a completely new region in 
which Am 2 is lower than the previous two ones, which we call very-low-LMA. If one chooses e' < 0, the preferred 
allowed region tends to higher values of Am 2 . 

Almost all our conclusions below are independent of specific sources of the non-standard neutrino interactions, that 
could be present in interactions with d-quarks, u-quarks or electrons. 

We have found that the main effects of the presence of the NSNI interactions are: 

• displacement of low-LMA region to lower (higher) values of Am 2 , for a positive (negative) value of e'. 



suppression of Earth regeneration at Am 2 ~ 10 5 eV 2 for positive values of e'. 



• Due to suppression of Earth regeneration, appearance of a new region of compatibility between solar and 
KamLAND data around Am 2 ~ 1CP 5 eV 2 , with no spectrum distortion for the low-energy SK and SNO data. 

• improvement of high-LMA fit quality for positive values of e' 

• a 1 kton-yr of KamLAND can make a strong statement about the existence of non-standard neutrino interactions. 
The striking signal of this NSNI would be the location of the prefered oscillation parameters in the very low or 
in the high LMA region. 

Note added: When we were finishing our paper, an article by Friedland, Lunardini and Pena-Garay 
I hep-ph/0402266) appeared, which discusses topics similar to the ones discussed in our paper, where we discuss 



not only the non-standard neutrino interaction induced by d-quarks case as well the u-quarks and electrons. Also 
we made a quantitative statement about the role of more statistics on KamLAND experiment, combined with the 
present solar neutrino data, to put more restrictive bounds on non-standard neutrino interactions. 
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Am 2 /4E (eV) 

FIG. 1: Survival probability of electron neutrinos. In this figure we can see the displacement in Am 2 /4E of the suppression pit 
associated with the transition between resonant and non-resonant regions. Also is possible to see the supression of regeneration 
effect for e' ~ 0.3. 
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FIG. 2: Allowed regions of oscillation parameters for 1 a, 90% C.L., 95% C.L., 99% C. L. and 3 a, using the constraints of 
solar+KamLAND data, with NSNI with d-quarks and where we minimized the x 2 with respect to the NSNI parameter e . The 
best fit point is marked by a star and happens for e' = 0. 
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FIG. 3: Future sensitivity of combined analysis of solar+KamLAND, assuming lkt-yr of exposure. The continuous line 
corresponds to the actual limit that is obtained with present KamLAND data. The dashed, dotted and long dashed curves 
refer, respectively, to the simulated parameters Am 2 and tan 2 6 lying at low LMA, very low LMA, and high LMA regions. 



